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Lobe Overflow as the Likely Cause of Pericenter Outburst in an 

SMBH Orbiter 

R.E. Wilson^’^ & E. J. Devinney, Jr.^ 

ABSTRACT 

A very large lobe overflow event is suggested to explain the 0.^4 brightening 
observed in K band at pericenter passage of the star known as S2 that orbits the 
Galaxy’s supermassive black hole (SMBH). Known observed properties of S2 that 
contribute to lobe filling are 1) the enormous mass ratio, Msmbh/Ms 2 , 2) S2’s 
fast rotation, and 3) S2’s large orbital eccentricity. Published estimates have given 
limiting lobe sizes of order 100 to 300 Rq but, with S2’s fast rotation taken into 
account, the computed lobe size is much smaller, being compatible with either 
a main sequence OB star or a stripped evolved star. An important evolutionary 
consideration that predicts very large pericenter overflows is envelope expansion 
following mass loss that is characteristic of highly evolved stars. Material removed 
by lobe overflow at pericenter is replenished by envelope expansion as an evolved 
star awaits its next pericenter passage. An observational signature of lobe overflow 
for upcoming pericenter passages would be appearance of emission lines as the 
ejected gas expands and becomes optically thin. 

Subject headings: Galaxy: center — stars: individual 


1. Background 


Prospects for understanding the evolut ion of the Galactic center region were boosted 
by discovery and quantitative investigation f Schodel. et ah 20021 ) of a star designated as S2 
that is bound to the Galaxy’s supermassive black hole (SMBH) in a 15.8 year or b it, wit h 
additional stars in the region under similar investigation - see also iGhez. et al.l fj2003al) . 
Many papers have since added observations and ideas on individual SMBH orbiters and 
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on the stati s tics an d c ollective properties o f stars i n the inner Galacti c center region, e.g. 


Ghez. et ah J2003b ch: 

Alexander & Morris (2 

0031: 

31um. et ah (200 

31: Eisenhauer. et ah 

(20051 

: Gillessen. et ah 

(20091: Davies & Kine 

(|20051 

: Gillessen. et ah 

(20131: Witzel. et ah 

(20141 

and an extensive 

review bv Genzeh et a 

L (I 2 OI 

0). 


Issues abound regarding t he evolutionary s t ates o f these stars and how t h ey came to be in 


the in n er Galactic center, e .g. iDavies fc Kind (1200511 : IZhang. et al.l (120131) : iGillessen. et ah 


(120091) : iGenzel. et ahl (1201011 . Stars arriving on nearly par abolic or even hyperb olic orbits 
could suffer major stripping on initial pericenter passage (IDavies fc Kingl 120051 1. with the 
lost material carrying away enough orbital energy to leave the remnant in an elliptical orbit, 
hlling its limiting lobe at pericenter. Alternatively, a star that has been trapped into a tight 
orbit while well detached from its lobe could later undergo evolutionary expansion and attain 
lobe hlling. Perhaps a lost binary companion may carry off the requisite orbital energy at 
hrst encounter with the SMBH and leave the remaining star bound. 


Eisenhauer. et ahl (120051 1 found most of the brighter inner orbiters to be BO to B9 main 


sequence stars, with S2 in the range 08 to B O, and with rotation velocities typical of B stars 
in the Galactic disk. IDavies fc Kingl (1200511 argued that they are actually tidally stripped 
remnants of AGB stars that now superhcially resemble main sequence stars and estimated 
S2’s mass at below a solar mass, specihcally about 0.8 Mq. All in all, mass estimates for S2 
that have been published or correspond to observed (main sequence) spectral types range 
from 0.8 to more than 20 Mq. Accordingly, we simply adopt 10 Mq for exploratory lobe size 
computations. 


2. Brightening of Star S2 Dne to Lobe Overflow 


S2 has continued as the most thoroughly discussed SMBH orbiter, largely due to its 
having been observed over more than a full o rbit, and having brigh t ened by about (l .^dO in 
K band, coincident with pericenter passage (IGillessen. et al.l 1200911 . IGillessen. et al.l offered 
seven ideas for explaining the brightening, however they then ruled out four of the ideas 
and argued against likelihood of the other three. Gonsequences of lobe overflow - a very 
common player in a variety of close binary issues - were not among the seven ideas. Lobe 
overflow might not be considered if one were thinking only of the huge orbital scale (of order 
1000 AU) compared to the size of a main sequence star, but it turns out that S2’s limiting 
lobe is actually similar in size to a main sequence star of, say, 10 Mq, as shown in (jU Lobe 
overflow is an attractive idea for S2 brightening - one could postulate that S2 exceeded its 
limiting lobe at the 2002 pericenter passage and ejected a strong puff of material that quickly 
expanded in vacuum so as to appear as a rather large cloud of brightly emitting gas. The 
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lobe size issue will now be addressed. 


3. Lobe Size Essentials 

To place the present work in context, we review the formal relations in estimates for (1) 
tidal radii and (2) limiting lobes. Although these two terms are quite distinct, both often go 
by the name ’Roche limit’ and some recent papers treat them as equivalent, thereby leading to 
considerable confusion and perhaps even wrong conclusions. Item 1, tidal radius, concerns the 
distance from a mass at which an idealized fluid mass (usually a small satellite) is disrupted 
when tidal stretching matches or exceeds the satellite’s cohesion due to self-gravity. The 
simple tidal radius concept considers test particles on the surface of a self-gravitating sphere 
of mass m and radius r at a distance d from an object of mass M. The test particles, located 
on opposite ends of a diameter of m on the line of centers, suffer a stretching force (surface 
to center) per unit particle mass due to external mass M of 2GMr/d^, assuming r to be very 
small. The effective compressional force per unit particle mass between surface and center 
that can result in a static conhguration is the object’s surface gravity, Gm/r^. Quantity d is 
the tidal radius and marks the distance from M at which a very small object is disrupted, 
so the hnal relation pertains to the test particles being arbitrarily close together. If r is not 
small compared to d, then the simple relation may still give the tidal radius approximately, 
with only gravitation considered, although the full relation is then more complicated. A flaw 
in this picture with regard to stretched objects of hnite size is that the satellite is presumed 
spherical, whereas tidally stretched stars lack front-to-back symmetry (i.e. have ’’teardrop” 
shapes). But more important for the case of star S2 and probably other SMBH orbiters is 
that rotation of m is not considered in the traditional tidal limit development. In summation, 
the tidal limit relation between d and r, d = r(2M/m)^/^, can be inverted with d set to the 
orbital separation to give roughly correct limiting size where the assumptions apply, namely 
for small, synchronously rotating satellites, but may be wrong by orders of magnitude for 
fast rotating stars such as S2. 

Item 2 is commonly called a Roche lobe, although it was not originated or even consid¬ 
ered by Roche, who did however consider a special case of the potential utilized today. As 
the idea is to specify a size limit set by the condition that material not be spilled from a star, 
the descriptive term ’limiting lobe’ serves well. The insight that spawned the concept came 


^In this context, ’synchronous rotation’ means that angular star rotation and mean orbital angular rota¬ 
tion are equal or, equivalently, that the star rotates once per orbit period. Note that there are other meanings 
of ’synchronous rotation’, each valid in its own context. 
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from lKuiperl fjl94lh . who realized that tidal force is an unnecessary complication with regard 
to the lobe size limit, as only one point, not two, need be considered, and only ordinary 
effective gravity at that point, not differences between two points, need be computed. The 
procedure is (step 1) to hnd the point along the line of centers where ordinary gravitational 
(not differential tide raising) forces due to M and m, along with local rotational force, add 
to zero. Material that is stationary in a frame that co-rotates with the star is not bound 
to the star at this special point, so an ejection nozzle forms. Asynchronous examples are 
now treated via a factor (see §1]) that alters the centrifugal t e rm without a ffecting the 
basic idea of locating the effective gravity null point flPlavedllQSSl: lLimberlll963h . A dehnite 
equipotential that dehnes the star surface passes through the special point of null effective 
gravity, so (step 2) numerically integrate the volume, Viobe, enclosed by that equipotential 
and thereby hnd the equivalent-sphere mean radius as 


Rr, 


= (34^ofee/47r)'/l 


( 1 ) 


Kuiper assumed synchronous rotation, which is the expected and observationally indicated 
case for very close binaries (due to tidal locking). 

Conditions that lead to small limiting lobes for SMBH orbiters are the enormous mass 
ratio, large orbital eccentricity, and - not previously emphasized in the literature ~ fast 
rotation. For hxed SMBH mas^, lobe size decreases with decreasing star mass, decreasing 
orbit size, increasing eccentricity, and increasing star rotation. S2 is known to be a fast 
rotator, while any kind of tidal locking would produce exceedingly slow rotation in view of 
the 15.8 year flGhez. et al.ll2003d) orbit period, Porb- The fastest locked rotations would be for 
locking to the pericenter orbital angular rate and give Pmf aroun d half a year for S2, whereas 
the measured lAo^sini is 220 ±40 km sec~^ flGhez. et al.ll2003d ). so there is no tidal locking 
of any kind. Whether the star’ s equator is aligned w' ith the orbit plane in not known but the 
orbital sini is 0.7040 ±0.0058 flGillessen. et al.ll2009[) so, under the assumption of alignment, 
Vrot ~ 312 ± 57 km sec~^. The corresponding angular rotation, assuming Peg = 5.0Pq, 
is ~ 7100 times the mean orbital angular rate, and rotational force becomes important in 
setting local effective gravity. It will be shown below that the problem of main sequence 
SMBH orbiters being too small to exceed their limiting lobes can disappear if they rotate at 
typical B star rates, as does S2. 


^We adopt Msmbh = 4.31 x 10® Mq (|Gillessen. et al.l 120091) and thus a mass ratio, Msmbh/Miq, of 
4.31 X 10® for a 10 Mq star. 
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3.1. The Eggleton Approximation 


Approxima tion formulas are often used to estimate lobe size, most commonly one by 


EggletonI (119831) that reproduces accurately computed mean lobe radii, based on the Kuiper 
logic, to better than 1 percent over the full range of mass ratio, from 0 to ocjfl. Although a 
rather small computer program can generate such lobe radii with negligible error, and some 
public binary star programs list lobe radii as incidental output, the Eggleton approximation 
has provided a one-line lobe calculation in many evolutionary programs where 1 percent 
accuracy may be sufficient. However note that the Eggleton formula is specifically for syn¬ 
chronous rotation and not meant for stars that rotate faster or slower than synchronously. 
It will give limiting lobe radii that are too large by orders of magnitude if applied to fast 
rotators such as star S2, and may be responsible for misleading conclusions where rotation 
rates are unknown. An algorithm that follows the Kuiper logic, enhanced to handle arbitrary 
rotation and eccentricity, is not difficult to program and avoids the approximations of fitted 
formulas. Most accurate limiting lobe computations now adopt Kuiper’s strategy, usually via 
one of the commonly used binary system light/velocity curve progr ams, although another 
option can be the collection of intricate approximation formulas by ISepinski. et al.l (120071) 
that account for asynchronism and eccentricity. 


4. Quantitative Estimate of S2’s Lobe Size 


Computation of a binary component’s limiting lobe geometry begins with solution for 
a point of null effective gravity along the line of star centers (x-axis), thereby locating the 
nozzle from which matter flows if the lobe is filled or slightly overfilled. The relevant equation 
for the S2 problem must account for orbital e ccentrici t y and the star’s rotation in addition 
to the gravity of both objects, as does eqn. 3 of Wilson ( 1979 ) for the derivative of potential 
in the x-direction, which is zero at the null point. That equation is 


'^Incidentally, we checked the 1 percent accuracy statement in lEsfgletonl (jl98,4) at the request of Prof. 
Eggleton, finding no discrepancies as large as 0.8 percent among 18 widely spread mass ratios, of which only 
two exceeded half a percent. Column 2 of Eggleton’s Table 1 (mean lobe radii from integrated volumes) was 
reproduced to all printed digits except for two differences of 1 in the last plac e. Column 6, the Eggleton 
approximation, was reproduced exactly. The checks were done with the WD ( Wilson fc DevinnevI 1971 1 
computer model. 

^The potential is a modified version according to the convention in Kopal ( 19591) . 
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dQ 

dx 


X 


+ 


q{D — x) 


(a;2 + 1/^ + {[D — xY + 


+ F^(l + q)x — qjD^ 


( 2 ) 


Rotation enters via a parameter F, the ratio of rotational angular velocity bJrot to mean 
(i.e. time-averaged) orbital angular velocity cuorb- Other input quantities are the component 
mass ratio (g = M2/M1), momentary separation of star centers {D), and x, y, z rectangular 
coordinates of a point at which dfl/dx, and subsequently O, are to be computed. Here S2 
is taken to be star 1 and the SMBH is object 2, so the mass ratio is a large number rather 
than its reciprocal. The unit for x, y, and z is D while D is in unit a, the semi-major 
axis of the relative orbit, in computations with equation [21 with F = 1 — e at periastron 
or pericenter. The location of the null point along the line of centers is found by setting 
y = z = Q and dfl/dx also to 0, setting the dimensionless angular rotation F and eccentricity 
e to values of interest, and then solving for x by numerical inversion (such as Newton- 
Raphson iteration). The potential at the null point then establishes the lobe surf ace’s 3- 
dimen sional form as an equipotential that includes the null point (see eqn. 1 of I Wilson 
(119791) for the generalized defining potential). The equipotential’s enclosed volume {V) can 
then be integrated numerically via the defining equation and a mean lobe radius found from 
eqn. [1] A final step computes equatorial rotation velocity, 14^, from angular velocity. That 
calculation is simplihed by the star being almost axially symmetric and its equator circular 
at these fast rotation rates, so there is no issue of where along the equator the result applies. 
Accordingly 


beg R'eq^orbFi (3) 

with length in km, time in seconds, and mean orbital angular velocity, ujorb = ‘^'^/Porbi in 
radians / sec. 

The binary star modeling and analysis program (WD program@) applied here has re- 
hnements that allow reliable operation in difficult circumstances. For example, its Newton- 
Raphson iterations (for inversion of equation[2]to hnd the effective gravity null point) evaluate 
several Taylor series terms beyond the usual first derivative term. This point is mentioned 
so that readers who may write their own inversion program to check our results are not 
disappointed by failed computations. A relatively simple inversion scheme can converge well 
for ordinary mass ratios but not for ultra-large mass ratios such as the 4.31 x 10® of the 


^The WD program’s most recent public version, with documentation and sample input files, can be 
downloaded from anonymous FTP site ftp.astro.ufl.edu. Go to sub-directory pub/wilson/lcdc2013. 
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present problem Also important for fractionally tiny lobes (large g, large F) is to begin iter¬ 
ations already close to the nnll point, so as to avoid an initial jnmp beyond the proper range 
between the star centers, from which recovery is difficult. Fortunately such a conhguration 
admits particularly good starting estimates of the null point’s location. To see this readily, 
write eqn. |2]as it applies along the line of centers at the null point. 

This form is a quintic equation in x, soluble only iteratively, but with x very small the 
second and fourth terms on the right side very nearly cancel so that the remaining terms 
(also replacing 1+q with the very large q) lead to a simple result. 


X 


^-2/3^-l/3_ 


(5) 


The approximation is reasonably accurate only for quite small x, although very accurate for 
SMBH orb iters and perhaps usefully accurate for M 2 /Mi of a few hundred or more. 


Inputs to the lobe size computation for S2 were e = 0.88 (iGillessen. et al.N2009f) and 
Msmbh/Ms 2 = 4.31 X 10® (mass ratio), along with a few well spaced F’s. One of the F’s 
is close to the nominal value of 7100 that goes with our rough estimate of Veg that assumed 
alignment of the equatorial and orbit planes in §21 The resulting mean lobe radius is 6.5Rq, 


which is larger than a 10 Mq main sequence sta r (about 3 on t 


le ZAMS to 5 Rq at the 


TAMS), although the spectral type estimate bv lEisenhauer. et al.l (120051) extends to 08, for 


which a main sequence radi us can exceed G.dRrv^. A stripped highly evolved star that resembles 
a main sequence star, as in iDavies &: Kind (120051) . remains a candidate. With either kind of 
star, the idea of lobe overflow at pericenter passage now becomes a real possibility. Table [H 
has mean lobe radij^ for four assumed angular rotation velocities (F’s) of the 10 Mq model 
orbiter to give a sense of how steeply lobe size depends on rotation rate. A check to see if 
the program gives the right order of lobe size is provided by calculation of the equatorial 
radius of a 10 Mq isolated star (no SMBH) that is marginally unbound at the equator 
while rotating at one of the table values, 307 km sec~^. If the magnitudes of rotational and 
gravitational force are then equated, the equatorial radius will be given by Reg = GM/V'^, 
which evaluates to 2O.2i?0 for a 10 Mq star. The corresponding mean radius will be smaller 
since Rpoie is smaller than Reg, so rotation alone produces a limiting size only about three 
times greater than do the combined effects of rotation and the SMBH gravity. The purely 


®Note that these are ’equivalent sphere’ radii, not distances to the effective gravity null point. 


















gravitational lobe radius for a slowly rotating star, with e = 0.88 and the present problem’s 
adopted masses, is ~ lOOi? 0 , so the effect of fast rotation on lobe size is not small. 


4.1. Why Such Large Scale Ejection? 


A remaining issue is why a huge puff would be ejected at pericenter passage. The ordi¬ 
nary context of lobe overflow is the synchronous-circular case that is commonly encountered 
in close binary systems, where gas leaks out quiescently and is usually difficult or impossible 
to detect photometrically. S2, being a very fast rotator, will not undergo the gentle process 
of the synchronous-circular case with its low ejection velocity. The supersynchronous case 
is very different, with an ejection velocity close to the star’s equatorial velocity, which is of 
order 300 k m sec~^ for our r nodel of S 2 . And why would a large amount of gas be ejected? 
Suppose the iDavies fc Kind (120051) proposal, that the close-in orbiters ar e tidall y strip ped 


highly evolved stars, is correct, and that S2 is typical. Well known (e.g. iPlaved (119681 )) is 


that radii of highly evolved (i.e. chemically stratihed) stars increase with loss of envelope 
matter, in contrast with shrinkage for unevolved and modestly evolved stars. S2 has 15.8 
years between pericenter passages to expand following each pass and could arrive at pericen¬ 
ter not just marginally hlling its lobe but substantially overhlling it. Although a quantitative 
estimate of the overhlling will require reasonably good estimates of S2’s internal structure 
that are not now in hand, the qualitative picture is that S2 may reach pericenter ready to 
send very fast moving gas through a large open nozzle, leading to a very large ejection event. 
One test of this idea, waiting for the next pericenter passage, is that emission lines should 
appear as the ejected gas expands and becomes optically thin. Naturally some or all of these 
expectations may be anticipated as other SMBH orbiters pass through their pericenters. 
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Table 1. Dependence of ’Equivalent Spherical Volume’ Lobe Radius on Equatorial 
Velocity for a 10 M©, e = 0.88 SMBH Orbiter Representing S2 


F (rotation parameter) 

Equatorial Velocity {km ■ sec 

’Equivalent Sphere’ Lobe Size {Rq 

1000 

44 

23.6 

3000 

131 

11.4 

7000 

307 

6.5 

10000 

438 

5.1 






